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Plate tectonics is largely responsible for material and heat circulation in Earth, but for unknown reasons it 
does not exist on Venus. The strength of planetary materials is a key control on plate tectonics because 
physical properties, such as temperature, pressure, stress, and chemical composition, result in strong 
rheological layering and convection in planetary interiors. Our deformation experiments show that crustal 
plagioclase is much weaker than mantle olivine at conditions corresponding to the Moho in Venus. 
Consequently, this strength contrast may produce a mechanical decoupling between the Venusian crust and 
interior mantle convection. One-dimensional numerical modeling using our experimental data confirms 
that this large strength contrast at the Moho impedes the surface motion of the Venusian crust and, as such, 
is an important factor in explaining the absence of plate tectonics on Venus. 

Venus is regarded as a twin planet to Earth because of the similar density, size, internal structure, and 
distance from the Sun of the two planets. However, the Magellan mission to Venus reported no evidence 
for a global pattern of plate tectonic margins and subduction on Venus, and that the Venusian surface has a 
nearly uniform age\ These observations imply that plate tectonics does not exist on Venus^ ^ although it is the 
primary control on material and heat circulation on Earth. Consequently, the absence of plate tectonics on Venus 
as compared with Earth might explain the different evolutionary histories of these otherwise similar rocky 
planets. 

A number of hypotheses have been proposed to explain the absence of plate tectonics on Venus, based on 
observational data and numerical modeling. Nimmo and McKenzie^ suggested that the lithosphere of Venus is 
unable to create faulting on the scale necessary to form plate boundaries due to the absence of water. In contrast, 
the large viscosity contrast between a planetary surface and mantle interior is a key factor for stagnant-lid 
convection, as when the planetary surface is stiff and sluggish as compared with the planetary interior, mantle 
convection can only occur beneath the lithosphere^ In these modeling studies, rheological structure has been 
shown to play an important role in the evolution of Venus. However, the rheological behavior of planetary 
interiors is sensitive to temperature, water, and chemical composition, which may all produce strong rheological 
layering^"^. Rheological decoupling may exist across the chemically distinct layers at the Moho where olivine-rich 
mantle is overlain by plagioclase-rich crust (Fig. 1). In this study, we carried out laboratory deformation experi- 
ments to constrain the nature of rheological layering across the Moho, and we discuss how this may explain the 
absence of plate tectonics on Venus. 

Rheological structure can be inferred from flow laws that represent the strength of solids, which are dependent 
on strain rate, temperature, and chemical composition^°'^\ Mackwell et al^ determined the flow law for diabase 
using rock deformation experiments. These results, when applied to rheological structures in Venus, suggest that 
a large strength contrast exists between its crust and mantle. In their model, rheological structure was modeled 
using the power-law creep formula: 



= A(7" exp 



RT 



(1) 



where e is the strain rate, A is a constant, q is the stress, n is the stress exponent, is the activation enthalpy, R is 
the gas constant, and Tis the temperature^^. This type of flow law is commonly applied to high-temperature creep, 
although the Peierls mechanism becomes dominant at low temperatures and high stresses^^. In the Peierls 
mechanism, the strain rate is exponentially proportional to applied stress as follows: 



— L-£{-(^)Ti 



(2) 



SCIENTIFIC REPORTS | 4 : 4403 | DOI: 1 0.1 038/srep04403 



1 



(a) Differential stress (MPa) (b) Differential stress (MPa) 

0 200 400 600 800 ^0 200 400 600 800 




Figure 1 | Strength profile models of Earth and Venus, (a) Strength profile of oceanic lithosphere on Earth. The flow law for the Peierls mechanism of 
olivine in the upper mantle^° and a relatively broad brittle-to-plastic transition (dashed line) were adopted from Kohlstedt et a?^. (b) Strength profile 
model for Venus predicted from our experimental results showing the large strength contrast between the lower crust and upper mantle. These 
models simulate deformation strength of the Venusian lithosphere at an age of 200 Ma, which is the estimated age of the surface around Artemis 
Chasma^^. In this model, the Moho depth was assumed to be the 20 km^^. Flow law for the Peierls mechanism^° and power-law creep"'^*^ were used for the 
mantle. Peierls mechanism for plagioclase, as determined in this study, was applied to the lower crust. The parameters of the flow laws used in this study 
are summarized in Table S3. 



where Op is the Peierls stress, and p and q are non-dimensional 
parameters that depend on the geometry of kinks Assuming a 
thermal gradient of 5-30 K/km and a Moho depth of 7-20 km^^, 
the Peierls mechanism is likely to be the dominant mechanism of 
deformation at the Moho in Venus. Therefore, the power-law creep 
formula determined by previous studies should not be applied to the 
Theological structure of the shallow part of Venus. In this study, we 
conducted two -phase deformation experiments to directly constrain 
the strength contrast between the crust and mantle at conditions 
corresponding to the Venusian Moho. 

Results 

Deformation experiments were carried out at a pressure of 2.0 GPa 
and temperatures of 600-1000°C using a sohd-medium deformation 
apparatus at Hiroshima University, Japan. The starting materials of 
annealed synthetic plagioclase (An66Ab34) and olivine (FogoFaio) 
were sandwiched separately between alumina pistons cut at 45° to 
the maximum compressional direction and oriented to represent the 
deformation geometry of simple shear (Fig. 82). The experiments 
were conducted under dry conditions at a constant rate from 6.5 
X 10"^ to 3.7 X 10"^ (s~^). Shear strain was measured from the 
rotation of a nickel strain marker that was initially orientated per- 
pendicular to the shear direction (Fig. 2). 

Samples recovered after the deformation experiments showed the 
development of a crystal-preferred orientation in plagioclase and 
undulatory extinction in olivine, which suggest that deformation of 
these minerals was largely controlled by inter- crystal deformation, 
such as dislocation- controlled creep. The shear strain determined 
from the rotation angle of the strain marker was up to y = 0.8 for 
olivine and y = 4.7 for plagioclase, and in nearly simple shear geo- 
metry (Table SI). The uncertainty on these shear strain measure- 
ments is ca. 10%-20% based on the shape of the strain marker, and 
the compressional strain was typically <20% of the shear strain 
(Table SI). The stress values should be identical between the two 
layers in our experimental assembly. Hence, the strain contrast pro- 
vides an approximation of the strength contrast at a given stress, 
temperature, and pressure. Figure 3 shows the strength contrast 



between plagioclase and olivine as a function of temperature. The 
strength of plagioclase was lower than that of olivine at any experi- 
mental temperature (filled square symbols) and strain rate (filled 
diamond symbol). These results are the opposite of those predicted 
from previous studies based on the power-law creep of plagioclase 
(dashed line), indicating that another deformation mechanism oper- 
ated under our experimental conditions. 

A large amount of extrapolation is needed when applying labor- 
atory data for plastic flow to geological processes, and this extrapola- 
tion must be made using an appropriate flow law for the deformation 
mechanism. Therefore, we determined the dominant deformation 
mechanism of olivine during our experiments with reference to a 
deformation mechanism map^^'^^"^\ which shows that the Peierls 
mechanism best explains our experimental results at temperatures 
below 1000°C (Fig. S5). Although the deformation mechanism of 
plagioclase is not well constrained, our experimental results at low 
temperatures clearly deviate from the power-law creep determined 
by previous studies^^'^^ in which plagioclase was shown to be stron- 
ger than olivine (Fig. 3). This indicates that the Peierls mechanism 
might also control plagioclase deformation under our experimental 
conditions. Tsenn and Carter^^ reported that the Peierls mechanism 
becomes dominant at low temperatures in materials with relatively 
strong chemical bonding such as silicates. As such, the strength 
contrast observed in our experiments is most likely caused by differ- 
ences in the Peierls stress of olivine and plagioclase. 

Discussion 

A strength profile was calculated by extrapolating frictional sliding 
and the viscous flow law of each material to high pressures and 
temperatures corresponding to a range of depths^. Figure 1 shows 
the strength profiles of Earth and Venus, and highlights that Earth's 
oceanic lithosphere has no strength contrast between crust and man- 
tle because deformation at the Moho is still in the brittle regime 
(Fig. la)^^. In this case, the crust and mantle are strongly and dynam- 
ically coupled and will move coherently together. In contrast, our 
experiments show a large strength contrast between the crust and 
mantle on Venus (Fig. lb), suggesting that the lower crust may act as 



SCIENTIFIC REPORTS | 4 : 4403 | DOI: 1 0.1 038/srep04403 



2 




10^ 



Strain marker 



(c) 



0.5 mm 





i 






Strain marker * 




^^^^^ 





Figure 2 | Experimental samples after shear deformation, (a) Olivine and 
plagioclase samples were sandwiched between alumina pistons, which were 
cut at an angle of 45° to the maximum compressional direction. The 
recovered sample was cut along the direction of compression, (b and c) 
The shear strain of the recovered samples was measured from the rotation 
of the strain marker that was initially perpendicular to the shear direction. 
Arrows indicate the sense of shear. 



a lubricant and therefore inhibit horizontal motion of the entire 
crust. In this case, the soft lower crustal material may be separated 
from the convecting mantle and never be able to be subducted into 
the stiffer mantle material. In support of this hypothesis, we note that 
Brown and Grimm^^ suggested that a landform called Artemis 
Chasma on Venus was produced by underthrusting and strike-slip 
deformation of the lithosphere, whereas the overlying crustal mate- 
rials were not subducted. The strength profile from our study indi- 
cates that this landform may have resulted from decoupling in the 
lower crust, which in effect forms a "crustal asthenosphere". On 
Earth, similar landforms are observed in continental collision zones, 
caused by complex deformation of the upper crust, including fold 
and thrust belts, and strike-slip faulting. The tectonic and rheological 
models in such zones emphasise the importance of decoupling 
between the crust and mantle, indicating that the overlying crust 
was not subducted^^. Shemenda^^ carried out analogue subduction 
experiments incorporating "crustal asthenosphere" for simulating 
continental collisions on Earth, and produced successive fold and 
thrust belts in the upper crust, in which the crust could not intrude 
into the mantle. We consider that a similar phenomenon could have 
occurred beneath Artemis Chasma in the absence of plate tectonics 
on Venus. The elastic thickness in Venus has been constrained to a 
range of 0-20 km in nearly half areas^^, although Artemis Chasma 
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Figure 3 | Strain contrast between plagioclase and olivine as a function of 
temperature. Square and diamond symbols denote the experimental 
results for strain rates of — 10"^ and —10"^ s~\ respectively. Dashed lines 
represent the Strain contrast calculated from Peierls mechanism for 
olivine^" and power-law creep of plagioclase^^'^^ at a strain rate of 10~^ s~\ 
The horizontal error-bars come from temperature gradient within 
samples. 

shows relatively thick elastic layer of 37 km^^. These estimates of 
elastic thickness maybe caused by the different depth of Moho where 
decoupling occurs and the overlain crust is nearly stagnant. 
Mechanical decoupling could cause the decrease of the elastic 
thickness^°'^\ 

To evaluate the influence of the strength contrast at the Moho on 
Venusian tectonics, we numerically analysed viscosity and moving 
velocity in the lithosphere using a one- dimensional model that simu- 
lates an embryonic subduction zone such as Artemis Chasma^^ 
(Supplementary Information; Fig. S3). In this analysis, simple shear 
deformation of the lithosphere beneath the thrust is calculated with 
fault friction at the lithospheric surface and the motion of the mantle 
lithosphere. The strength profiles shown in Figure 1 were used for the 
rheology models of the Venusian lithosphere. The Moho depth in 
our model was assumed to be 20 km based on gravity and topo- 
graphy data^^'^^ (Fig. lb), although the numerical calculations were 
also conducted at various crustal thicknesses (20, 45, and 70 km) and 
produced similar results in each case. Various models with the crust- 
mantle viscosity contrast ranging from 1 to 10^ were tested. Under 
these conditions, the resultant viscosity and moving velocity depth 
profiles are shown in Figure 4. The velocity profile demonstrates that 
crustal motion becomes more with increasing viscosity contrast at 
the Moho (Fig. 4). Our modeling suggests that a large viscosity dif- 
ference between the crust and mantle results in a more or less neg- 
ligible motion of the crustal surface, even though the plate faulting 
crosses the entire brittle deformation range and the underlying man- 
tle moves at a velocity of 10 cm/yr (Fig. 4). The crustal thickness in 
Venus is not well constrained; however, our models with different 
crustal thicknesses results in a similar moving velocity profile in each 
case (Fig. S4). Although the experimentally determined strain con- 
trast at the Moho is relatively small (—one order of magnitude), once 
decoupling between the crust and mantle initiates, the motion of the 
upper crust becomes sluggish and the stress in the crust decrease; 
consequently, a positive feedback could accelerate the decoupling 
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Figure 4 | Viscosity (left) and moving velocity (right) profiles calculated 
using a one-dimensional deformation model for the Venusian 
lithosphere. The numerical modeling results based on our experimental 
data show that the planetary surface becomes immobile when there is a 
large strength contrast at the Moho. (a) Model with no strength contrast at 
the Moho. (b and c) Model with a viscosity contrast of 10^ and 10^ Pa s at 
the Moho, respectively. The internal velocities obtained for each model are 
shown assuming that the crustal thickness is 20 km and the bottom 
(mantle) velocity is 10 cm/yr. In case c, the velocity is truncated at zero as 
the negative value of the crustal moving velocity means that the crust 
becomes immobile with a smaller fault friction than that used in this study. 

between the crust and mantle. This results in a convection style with a 
stagnant crust that is apparently similar to stagnant-lid type convec- 
tion^^. Although Aitta^^ suggested that the small Ra number of Venus 
played a key role in shaping the different tectonics of Earth and 
Venus, a viscosity contrast at the Moho facilitated the rheological 
decoupling that contributed to the absence of plate tectonics on 
Venus. 

In the past, Venus may have had a higher potential temperature 
than at present. If the potential temperature was 150-200°C higher 
than at present, which is assumed to be similar to that on Earth 
during the Archean^^, the crustal portion would have been weaker 
than the mantle (Fig. 3). Therefore, decoupling could have occurred 
at the Moho during Venus' past, suggesting stagnant-lid type con- 
vection. However, if the crustal thickness was substantially thinner in 
the past (< ~7 km), deformation in the lower crust would have been 
controlled by frictional sliding, meaning that a strength contrast at 
the Moho would have been unlikely to exist. In the later case, the 
crust might have been coupled with the mantle convection and sub- 
ducted along with the mantle lithosphere. 

We conclude that rheological structure and, in particular, the 
strength contrast at the Moho, might be a critical factor that has 
shaped the different tectonic evolutions of Earth and Venus. Note 



that an application of our laboratory results to the Venusian tectonics 
requires a large extrapolation; consequently, an understanding of the 
active deformation mechanism is key to understanding the tectonic 
evolution of terrestrial planets. 

Methods 

Deformation experiments were carried out at a pressure of 2.0 GPa and temperatures 
of 600-1000°C using a solid-medium deformation apparatus at Hiroshima 
University, Japan. The starting materials of annealed synthetic plagioclase 
(An66Ab34) and olivine (FogoFaio) were sandwiched separately between alumina 
pistons cut at 45° to the maximum compressional direction and oriented to represent 
the deformation geometry of simple shear (Fig. S2). The experiments were conducted 
under dry conditions, and the water contents of the starting materials were verified by 
Fourier-transform infrared spectroscopy as being ca. <30 ppm H/Si for both 
minerals (Fig. SI). Tiny grooves were made at the interface between the sample and 
piston to prevent slip during deformation. A nickel jacket surrounded these materials, 
and the oxygen fugacity was buffered by the Ni-NiO reaction. The pressure was first 
raised to 0.1 GPa, and then the temperature was increased at a rate of ~15°C/min to 
400°C. The pressure was then increased to 2.0 GPa and the temperature was 
increased to the desired value. Temperature was monitored by two Pt/Rh thermo- 
couples placed close to the sample. When pressure and temperature had attained the 
desired values, a piston was advanced at a constant rate from 6.5 X 10"^ to 3.7 X 10"^ 
(s" ^). The samples were rapidly quenched by switching off the thermocontroller when 
the deformation experiments were finished, and the pressure was then reduced. Shear 
strain was measured from the rotation of a nickel strain marker that was initially 
orientated perpendicular to the shear direction (Fig. 2). 
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